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ABSTRACT 

We study the spinning down time scale of rapidly rotating white dwarfs (WDs) in 
the frame of the core-degenerate (CD) scenario for type la supernovae (SNe la). In 
the CD scenario the Chandrasekhar or super-Chandrasekhar mass WD is formed at 
the termination of the common envelope phase or during the planetary nebula phase, 
from a merger of a WD companion with the hot core of a massive asymptotic giant 
branch star. In the CD scenario the rapidly rotating WD is formed shortly after the 
stellar formation episode, and the delay from stellar formation to explosion is basically 
determined by the spin-down time of the rapidly rotating merger remnant. We find 
that gravitational radiation is inefficient in spinning down WDs, while the magneto- 
dipole radiation torque can lead to delay times that are required to explain SNe la. To 
explain the delay-time-distribution of SNe la the merger remnants distribution should 
be dN /dlog(B s'm 5) ~ constant, for 10 6 G < Bsin5 < 10 8 G where B is the dipole 
magnetic value and 5 the angle between the magnetic dipole axis and rotation axis. 



1. INTODUCTION 



Type-la supernovae (SNe la) are widely thought to be the thermonuclear detonations of carbon - 
oxygen white dwarfs (WDs) whose masses are near the Chandrasekhar limit (jHovle &: Fowlerl ll960), 
The route to a chieve this mass and the condit ions for explosion are in dispute (see reviews by, e.g., 
Livid l200ll and iHillebrandt &: Niemeverl 120001) . Two main routes are discussed in the literature. In 
the s ingle degenerate (SD) scenario (e.g., Whelan k, ibenfl973 : Nomoto 1982 ; Han k, Podsiadlowski 
2004) a WD grows in mass throu gh accretion from a non-degenerate st ellar companion. In th e 



doub le degenerate (DD) scenario (jWebbink 



1984 



Iben fc Tutukov 



1984; see 



van Kerkwiik et al 



2010 for a recent paper on sub-Chandrasekhar mass remnants) two WDs merge after losing energy 
and angular momentum through the radiation of gravitational waves. Observations and theoretical 
studies cannot teach us yet whethe r both models for SNe la can work, only one of them, or none 
(see recent summary bv lMaoal2010l ). 



The observations and theoretical results that favor and disfavor eac h scenar i o for the progeni- 
tors of SNe la are summarized by, e.g.. iLiviol (|200ll ) and more recently bv iHowelll (|201ll ). The main 
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problem for the DD s cenario (see review bvlHoweilll201lr ) is that in many cases an off-center carbon 
ignition occurs (e.g., ISaio k, Nomotd 12004 ) leading to accretion indu ced collapse (AIC) where a 
neutron star (NS) is formed rather than a SNe la. lYoon et al.l (|2007l ) raised the possibility that 
in a merger process where the more massive WD is hot, off-axis ignition of carbon is less likely to 
occur. The reason is that a hot WD is larger, such that its potential well is shallower and the peak 
temperature of the destructed WD (the lighter WD) accreted material is lower. Hence, in such a 
case the supercritical-mass remnant is more likely to ignite carbon in the center at a later time, 
leading to a SN la. Namely, the merger remnant becomes a rapidly rotating massive WD, that can 
collapse only after it loses sufficient angular momentum. 



Motivated by the possibility raised by lYoon et al.l (|2007l ). Kashi fe Sokerl (|201ll ) suggested 
recently that some SNe la might originate from a double degenerate merger that occurs at the end 
of the common envelope (CE) phase, rather than a long time after the CE phase as in the canonical 
DD scenario. Practically, the merger is of a companion WD with the core of an AGB sta r. The 
merger of a WD with the core of an AGB star was proposed by ISparks &: Stecherl ()1974l ) . They 
consider the outcome to be a SN explosi on that forms a NS, e.g., type II SN explosion. More 
relevant to this paper is the suggestion of iLivio Sz Riesa (|2003l ) that the merger of the WD with 
the AGB core leads to a SN la that occ urs at the end of the CE phase or shortly after. In that 
case hydrogen lines will be detected. As ILivio fc Riessl (|2003l ) aim was to explain the presence of 
hydrogen i n a SN la, they did not consider a delayed explosion. In this paper we discuss delayed 
explosion. Livio fc Riessl (20031) noted also that for such a a merger to occur the AGB star should 



be massive (also lKashi &: Sokerl 120 111 ). Tout et al. (2008) consider a merger of a WD with a core 
of an AGB star to explain t he fo rmation of massive rotating WDs with strong magnetic fields. 
Wickramasinghe fc Ferrarid (|2000l ) commented that such WDs might be more likely to form SNe 
la. 

Either the companion WD is destructed and forms a disk around the hot core, or the more 
massive core is destructed (because it is hot and large) and forms an accretion disk around the WD. 
In the later case the lower mass WD has a shallower gravita tional potential , hence the temperature 
in the accretion disk does not reach ignition temperature. iKashi Sokerl (|201ll ) termed this the 
core-degenerate (CD) scenario. In the CD scenario the merger occurs in a CE with a massive 
AGB star, or shortly (< 10 5 yr) after the CE phase, i.e. during the planetary nebula phase . 
Because of its rapid rotation the super-Chandrasekhar WD does not e xplode (e.g.. lAnandlll965l : 
Ostriker &; Bodenheimerlll968l ; lUenishi et al.ll2003l ; lYoon &: Langerll2005l ) . It will explode long after 
the merger, only after it has spun down to allow explosion. The CD scenario, including findings 
from this paper, is summarized schematically in Figure [TJ This scenario is completely different from 



the m odel of an explosion inside a CE wit h a low ma s s evol ved star as proposed by lHachisu et al 
(| 19891 ) that was subsequently criticized by lApplegatd (jl99ll ) . 



There are three key ingredients in the CD scenario, in addition to the common condition that 
the remnant mass be > 1.4M Q . (1) The hot core is more massive than the companion cold WD. In 
a future paper we will conduct a population synthesis to find the expected number of such systems. 
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The Core-Degenerate Scenario for SNe la 
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Fig. 1. — A schematic summary of the core-degenerate (CD) scenario for SNe la (jKashi 8z Soker 
201 ll ). This paper deals with the spin-down process. 
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(2) The merger should occur while the core is still large, hence hot. This lim its the merger to 
occur within ~ 10 5 yr after the common envelope phase. iKashi &: Sokerl ( 201ll ) showed that this 



condition can be met when the AGB star is massive, and some of the ejected CE gas falls back. (3) 
The delay between merger and explosion should be up to ~ 10 10 yr if this scenario is to account 
for SNe la in old-stellar populations, in addition to SNe la in young stellar populations. This is 
the subject of this paper. 

In section [2] we examine the spinning-down time by gravitational waves. Contrary to previous 
claims, we find that this time is very long, and might not be relevant to spin down WDs, despite 
being very efficient in spinning down neutron stars (NS). In section [3] we show that it is quite 
plausible that magneto-dipole radiation torque can lead to a long delay of the explosion, up to 
~ 10 10 yr. In section H] we constrain the properties of the core- WD merger product if the CD 
scenario is to account for a large fraction of SNe la. Our discussion and summary are in section [5j 



THE LIMITED ROLE OF GRAVITATIONAL RADIATION 



As done by many authors (e.g., 


Andersson et al. 


1999, 


Yoon & Laneer 


2004 




Yoon & Laneer 


2005 




Yoon et al. 


2007 




Justham 


2011: 


Di Stefano et al. 2011 


) we consider the emission of grav- 



itational waves due to r-modes in rotating WD progenitors of SNe la, but contr ary to most of 
these and other works, we will find this mechanism to be inefficient. We note that Piersanti et al 



(|2003al lb|) considered gravitational waves from WDs that have rotational to binding energy ratio 
of T/W = 0. 14, and therefore are highly deformed. The WDs we consider cannot reach values of 



T/W = 0.14 (jYoon fc Laneer! 12003 ). 



We follow the derivation given bv lLindblom et al.l (j 19981 ). where the assumptions and approx- 
imation can be found. We will take the mode I = m = 2, which is the most significant one. The 
velocity perturbation of the r mode is 



51? = aRfl ffie^ 
where Y^ is the magnetic type vector spherical harmonic, the frequency is given by 



l + l 



-n, 



(i) 



(2) 



fi and R are the stellar an gular velocity an d radiu s, respectively, and a is the amplitude of the mode. 
From equations 14-17 in iLindblom et al.l (|1998l ) we derive the power, Wgw> of the gravitational 
radiation (for / = m = 2, and written in a form that can be traced to their formulae) 

Wow = ( I) 7 I* P (rVdr\ 2 iT W. (3) 



We define the structural constant f3 by 



/ p{r)r 6 dr = @MR 4 

Jo 287T 



(4) 
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where (3 = 1 for a constant density sphere of mass M and radius R. We also define the ratio of 
the angular velocity to the break-up angular velocity (the Keplerian angular velocity on the stellar 
equator f^Kep) 

= - = - ■ (5) 



^Kep {GM/R?)W 
With these definitions the gravitational radiation power becomes 

2 13 q5 



Wqw 



7 2 x 5 2 x 3 6 tt c 7 



M 6 R- 6 Q 8 f3 2 a 2 . 



The moment of inertia is parameterized with the structural constant /?/, 

J = OAfaMR 2 , 



(6) 



(7) 



where for a constant density sphe re j3i = 1. We integrate for the structure of a massive WD model 
taken from lYoon Langerl (|2005l ) . and find f3 ~ 0.12, and f5i ~ 0.27. Assuming that all the power 
of the gravitational waves comes from the rotational kinetic energy, the spin-down timescale of 
the WD is given by integrating the equation d(I£l 2 /2)/dt = Wgwi where a solid body rotation is 
assumed here. This gives the timescale for spinning down from an initial fast rotation (Iq down to 
a critical (scaled) angular velocity Cl c , 

/ ~ \ -6" 

/ n 



TGW 



ml 



6W GW (n c ) 



1 Uc 



10 



11 



10 



M 
1.5M 




a 



10 



-2 




-6 



(8) 



yr. 



In spinning-down rapidly rotating NSs, that have a much smaller radius than WDs, gravita- 
tional radiation emitted by the r-modes is very efficient. For the gravitational radiation to play any 
role in spinning-down WDs the amplitude of the r-modes must become a > 0.1. This is quite large. 
Here we raise the possibility that magnetic fields limit the growth of th e r-mode to a degr e e tha t 
makes the gravitational radiation insignificant in spinning down WDs. ICuofano fc Dragol (|2010i ) 
find that as the magnetic fields in NSs become stronger, the lower limit on the rotation frequency 
for the development of r-modes instability becomes lar ger. In other cases t he ma gnetic field limits 
the amplitude of the r-modes. In the calculations of ICuofano &; Dragol (|2010i ) an equipartition 
between the magnetic pressure and the mode energy density in the outer region is achieved where 
the density is about one per cent of the average density p ~ 0.01p av . Using the same approximate 
equipartition, (l/2)p(aQ ro tR) 2 ~ B 2 /8tt, we obtain 

-2 




M 



1.5 M Q 



G. 



(9) 



This is a huge magnetic field for a WD. Before reaching this value the magnetic field will dissipate, 
and is likely to prevent the growth of the r-mode amplitude to a value of a ~ 0.01. Our conclusion 
is that most likely gravitational waves from r-modes are not significant in spinning down WDs. 
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VERY LONG DELAYED EXPLOSION BY MAGNETIC BREAKING 



A misalignment between the magnetic axis and the rotation axis leads t o a spin-down process 



due to magneto-dipole radiatio n torque that is commonly used for pulsars fe.g..lContopoulos &: Spitkovskv 



2006 



Benacquista et al.l 120031 ) . This mechanism was used for WDs by IBenacquista et al.l (|2003l ) 



with the following expression for the spin-down time from an initial fast rotation Qq down to a 
critical angular velocity 0, c , 



Ic 3 




-2 



R 



4000 km 



sin 5 

~b~T 



(sin 5) 



0.3 



-2 



10 s 



B 



10 8 G 



-2 



a. 



-2 



0.7ft 



Kep 




(10) 



yr, 



where 5 is the angle between the magnetic axis and the rotation axis, and 13 j is defined in equation 
([7]). The scaling of ft c at which the WD expl odes is based on the results of lYoon fc Langerl ([2005) 
for WD in the lower range (1.4 — 1.5M ; also lOstriker &: Bodenheimerlll968l ). The angular velocity 
distribution is not considered here in detail; only in calculating the total angular momentum a solid 
body rotation is assumed. In a future paper this assumption will be relaxed. H owever, we do note 
that magnetic fields might enforce a solid body rotation in the WD (|Piroll2008l ). 



Stabilizing rapidly rot ating super-Chandra sekhar WDs is a delicate matter (e.g- jYoon Sz Langer 



2004. IChen &: Lill2009l . and lHachisu et al.ll201ll in the single degenerate scenario). The strong mag- 
netic fields required in the present model most likely will enforce a rigid rotation within a short 
time scal e due the WD being a perfect conductor. The critical mass of rigidly rotating WDs is 
1.48Mq (jYoon fc Langerl 120041 and references therein). This implies that in the case of a strong 
dipole filed that exists in most of the their interior, WDs more massive than 1.48M will explode 
in a relatively short time. The exact time will be estimated in a future study. The similarity 
of most SN la suggests that their progenitors indeed come from a narrow mass range. This is 
~ 1.4 - 1.48M Q in the CD scenario. 

Equations (|SJ), ([9|), and (fit)]) bring us to consider the following scenario. Massive WDs, A%d > 
1.5Mq that can explode (or collapse) with higher angular momentum, do so on a typical time scale 
of < 10 9 yr; most within a much sorter time. WDs in the lower mass domain (Mwd — 1.4— 1.48M ) 
that explode with lower angular momentum and that on average have weaker magnetic fields might 
survive for a time of > 10 9 yr from their formation in a merger process. The exact time of explosion 
depends strongly on the magnetic field of the WD and the inclination angle. We cannot constrain 
the value of these parameters and the evolution of the magnetic field during the spin-down process, 
as, for example, r-modes might amplify the magnetic field. However, in the next section we crudely 
constrain the distribution of the relevant parameters. 



-7- 



DELAY TIME DISTRIBUTION 



In the previous two sections we have studied the evolution of a massive rapidly rotating WD 
that is a merger product of two WDs, or an AGB core with a WD, termed CD scenario for SNe la. 
The rapidly rotating WD does not collapse because of the centrifugal forces. On ce it spins down 
below the critical angular velocity f2 c , it becomes unstable and explodes as a SN la (|Yoon &: Langer 
2005). We now check the possibility that most of the SN la are formed by the CD scenario, and 
use observations to constrain the properties of these SN la progenitors . 



Equation (jlOj) for the spin-down time depends on 4 physical variables. The radius R and the 
critical angular velocity Cl c will not differ much from one WD to another, as we are mainly interested 
in WDs of masses ~ 1.5M that are rapidly rotating. The magnetic field and the inclination angle 
can vary by more than two orders of magnitude between different WDs. For that, we write the 
spin-down time from equation (|10p in the form 



t b =F(n c ,R)r ] - 2 
where we defined what we term the magnetic-dipole parameter 

r] = B sin 5. 



(11) 



(12) 



We now try to use observations to estimate the distribution of WDs that are the remnants of 
core-degenerate mergers with respect to rj, under the assumption that most (all) SNe l a are formed 



throu gh this channel. Maoz et al. (2011; see also lGraur et al.ll201ll . iMaoz et al.ll2010l . iRuiter et al 
2010) summarize the SN la rate versus time — the delay time distribution (DTD) — as found from 
observations 



dN 
~dt 



oc t e 



(13) 



with e ~ — 1, and where t is the time since the star formation event that formed the binary system 
that later became the SN progenitor. The last equation can be written as 



dN 
~~dt 



dN drj 
drj dt 



From equation (jlip with tb — > t, we have 



dn 8 



that when substituted into equation (fT4"|) gives 



dN 3 , , 9 C _1 

— — oc i2+ e oc 7? 6 ze ~ n \ 
dr\ 



(14) 



(15) 



(16) 



We recall that this result holds only for the massive WDs that were formed from the merger of 
two lower mass WDs during or shortly after a common envelope phase under the assumption that 
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most SNe la are formed via the CD channel. Equation (|16p can be also written explicitly as 

dN 



constant 



for 



10° G < Bsin5 < 10 8 G, 



(17) 



dlog(B sm5) 

where the range o f rj = B sin 5 comes from the range of 10 7 < t < 10 10 yr of the observed SNe rate 



( Maoz et al 



2011 



and equation (|10p . We note that magnetic fi elds in massive WDs might be 
> 10% of their initial value at age of 10 10 yr ( Muslimov et al J 1 19951 ). If true, then the decay of the 
magnetic field should be taken into account for long delay times, but it does not pose a problem 
to the model. 



Let us compare this results with observations. IWickramasinghe &: Ferrarid (|2000l ) review the 
subject of magnetic WDs. Magnetic WDs have magnetic fields in the range 3 x 10 4 - 10 9 G. This 
overlaps with the value assumed here. Although many of them are called rapid rotators because 
they are much faster rotators than most WDs, they are much slower than the spin rate required 
here. It is quite possible that super-Chandrasekhar magnetic WDs can be stabilized even if the 
surface s pinning rate is lower t han that assumed here, if a differential rotation is maintained inside 
the WD (jYoon fc Langerlboosh . 



Another encouragin g observational result is that ma gnetic WDs tend to be more massive 
that non- magnetic WDs (jWickramasinghe &: Ferrario 20 00) . Indeed, o ut of the 16 magnetic WDs 
with well determined mass reported by IWickramasinghe Ferrarid (|2000l ). two have masses of 
M > 1.3M©. This fraction is further discussed in the next section. 



5. DISCUSSION AND SUMMARY 



In this paper we examined the delay time from merger to explosion of the core-degenerate 
(CD) scenario for SNe la. A key ingredient in this scenario is that the merger occurs while the 
more massive degenerate component is the core of an AGB star or a planetary nebulae, e.g., at the 
end of the common envel ope phase or short l y afte r Earlier studies of the merger of a WD with the 
core o f an AGB star are ISparks &: Stecherl (|1974l ) as a way to make a neutron s tar, iLivio &: Riess 
(|2003l ) as an evolutio nary route for SNe la with hydrogen lines, iTout et al.l (|2008l ) as a way to make 
magnetic WDs, and iKashi Sokerl (|201ll ) in relation to SNe la with long delay times as studied 
here. 

Such a merger where the mor e massive WD is s till hot, hence larger, is more likely to avoid an 
early off-center ignition of carbon (|Yoon et al.ll2007l ): an early off-center ignition leads to a NeMgO 
WD instead of to a SN la. The problem of off-center carbon ignition during merger of cold WDs 



is the most severe problem of the double degenerate (DD) scenario for SNe la (|Howellll201ll ), The 



CD scenario avoids this problem. The remnant of the core-degenerate merger is a rapidly rotating 
WD. If its mass is super- critical it will explode after spinning-down to a critical angular velocity 
f? c (Yoon &: Langerl 2005 ). and not necessarily in a very short time delay as in the scenario of 
Livio fc Riesl (fcoosh where one is expect to observe hydrogen lines. 
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In section [2] we examined the spinning-down time by the commonly used gravitational waves 
mechanism. Although this mechanism is know to be very efficient in spinning down neutron stars, 
we found that this mechanism is not efficient in spinning down WDs (eq. EJ) . It can be efficient only 
if the amplitude of the r- modes becomes a > 0.1. This is a very large amplitude and the modes 
are likely to be damped by magnetic fields before reaching this value (eq. [9]). On the other hand, 
in section [3] we estimated the spinning-down time due to magneto-dipole radiation torque (eq. [TO]) , 
and found it to operate over the required time scale of ~ 10 6 — 10 10 yr for plausible parameters of 
the massive rapidly rotating WD remnant of the core-degenerate merger. 

If the CD scenario for SNe la is to explain most (all) SNe la, it should account for the delay time 
distribution (DTD), i.e., the S Ne la rate v ersus time since their progenitor formation. Using the 
observations as summarized bv lMaozl (|2010l ). we derive a crude distribution of the dipole parameter 
rj = B sin 5, where B is the magnetic field and 5 the inclination angle between the magnetic dipole 
axis and the rotation axis. This is given in equation (|17p . 

Before the CD scenario can be accepted as one of the channels to form SNe la more detailed 
calculations should be done. We end with a summary of these calculations, and with some possible 
promising properties of this scenario. 



1. The claim of lYoon et al.l (|2007l ) that when the massive WD in a merger is hot (in the CD 
scenario it is the core of an AGB star or of a planetary nebula), the off-center ignition of 
carbon is likely to be avoided, should be studied in more detail. This study should include 
the process where the more massive core is destructed (because of its larger radius), and the 
lower mass WD becomes the core of the remnant. 

2. A population synthesis study is required to determine the number of core- WD mergers that 
lead to Chandrasekhar and super-Chandr asekhar mass r e mnan ts. This is a subject of a 
future paper. The first step was carried bv lKashi &: Sokerl (|201ll ) who argued that the core- 
WD merger can be a common outcome of the common envelope phase with massive AGB 
stars, ~ 5 — 8M . Tout et al. (2008) estimated that about three times as many common 
envelope events lead to a merged core as to a cataclysmic variable. This is a high fraction. 
Wickramasinghe &: Ferrarid (|2000l ) estimated that ~ 5% of all WDs are magnetic WDs. Even 



if a small fraction of these merger products lead to SNe la, the CD scenario can account for a 
large faction of SNe la. The exact numbers will be derived in a forthcoming paper. The large 
fraction of massive stars that are in binary systems gives us hope that indeed the number of 
such systems might be adequate to explain most (all) SNe la. We note here that SNe la from 



WDs masses much above the Chandrasekhar mass, as observed in some cases, e.g. iHowell 
(|200ll ). are accounted for in the CD scenario. 



3. A detailed study of the spin-down time of the merger remnant due to the magneto-dipole 
radiation torque is required. Here we showed that this can be a promising mechanism. The 
evolution of the magnetic field with time should be examined as well. 
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4. It is expected that less massive remnant (but still super-Chandrasekhar) will take longer time 
on average to spin down and explode. T his might acc ount for the finding that SNe la in 



older populations are less luminous (e.g., iHowelll l200ll ). This should be studied after the 



magneto-dipole radiation torque process is better explored. 



5. Rapi dly rotating WDs have shallow density profile near the equatorial plane (|Yoon &: Langer 



20051 ). Shallow density pr ofile favors incomplete silicon burning that is required by observa- 



tions (jPakmor et al.ll2010l ). This process deserved to be studied. 



6. Eventually, rapidly rotating massive WDs should be found. Claims for massive rapidly 
rotating and s t rongly mag netized WDs have already been m ade, e.g. see discussion by 
Malheiro et al.1 (|201ll ) and IWickramasinghe fc Ferrarid (|2000l ). We can crudely estimate 



the expected number of such WDs if the CD scenario account for most of SNe la, in the 
following way (Maoz, D., private communic ation 2011). T he SN la rate per unit stellar mass 
in Sbc galaxies is ~ 1.5 x 10~ 1 3 yr^M" 1 ()Li et al.ll201ll ). From the loc al stellar density o f 
O.O85M pc" 3 (|McMillaJ201lT ) and the local WD density of 0.0046 pc" 3 (jllarris et alJbood ). 



a number of stellar mass per WD of 18.5M0WD" 1 in the solar neighborhood is obtained. 
If these WDs are to explode within ~ 10 10 yr with a time time delay distribution of t~ l in 
the range 10 7 — 10 10 yr, then their average age is ~ 1.5 x 10 9 yr, and the fraction of SN la 
progenitors should be ~ 0.4% of all WDs. Namely, in a sample of 10 4 WDs in the solar neigh- 
borhood the n umber of rapidly rotating massive WDs should be ~ 40. As ~ 5% of all WDs 
are magnetic (jWickramasinghe fc Ferrarid 12000 ) . it is expected that ~ 10% of the magnetic 
WDs will be massive enough to become SNe la. This is a predicti on of this model. Indeed, out 
of the 16 magnetic WDs with well determined mass reported by IWickramasinghe fc Ferrario 
(2000), two have masses close to the Chandrasekhar limit, M > 1.3M@. The study of large 
samples of WDs is the subject of a future paper. 
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